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Abstract 

The existing formulations of the Schrodinger interpolating dynamics, which 
is constrained by the prescribed input-output statistics data, utilize strictly 
positive Feynman-Kac kernels. This implies that the related Markov diffusion 
processes admit vanishing probability densities only at the boundaries of the 
spatial volume confining the process. 

We extend the framework to encompass singular potentials and associ- 
ated nonnegative Feynman-Kac-type kernels. It allows to deal with general 
nonnegative solutions of the Schrodinger boundary data problem. The result- 
ing stochastic processes are capable of both developing and destroying nodes 
(zeros) of probability densities in the course of their evolution. 
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1 From positive to nonnegative solutions of para- 
bolic equations 

We continue an investigation of the celebrated Schrodinger's boundary data problem, 
_ Qi °f reconstructing the "most likely" evolution which interpolates between the 
prescribed input-output statistics data (usually analyzed in terms of nonvanishing 
probability densities) in a fixed finite-time interval, interpreted as a duration time 
of the process. 

In the present paper we focus our attention again on stochastic Markov processes 
of diffusion-type (see Ref. [JTTJ for a jump process alternative), which are associated 
with the temporally adjoint pair of parabolic partial differential equations: 

d t u(x, t) = Au(x, t) — c(x, t)u(x, t) (1) 

d t v(x, t) = —Av(x, t) + c(x, t)v(x, t) 

Here, c(x, t) is a real function (left unspecified at the moment) and the solutions 
u(x, t), v(x, t) are sought for in the time interval [0, T] under the boundary conditions 
set at the time-interval borders: 

Po(x) — u(x, 0)v(x, 0) (2) 

Pt(x) = u(x,T)v(x,T) 



p (x)dx = p (A) , / p T (x)dx = p T (B) 

A JB 

We assume that p is a probability measure with the density p(x), and A, B stand 
for arbitrary Borel sets in the event space. In the above, suitable units were chosen 
to eliminate inessential in the present context (dimensional) parameters, and the 
process is supposed to live in/on R 1 . 



As emphasized in the previous publications, [pf-|nfl, the key ingredient of the 
formalism is to specify the function c(x, t) such that exp[— J * H(r)dr] can be viewed 
as a strongly continuous semigroup operator with the generator H(t) = —A + c(t), 
associated with the familiar |12| Feynman-Kac kernel: 



(f,exp[- H(r)dr}g) = J dy J dxf(y)k(y,0,x,t)g(x) = (3) 



f(u(Q))g(u(t))exp[- / c(u(t), T)dr\dp {uj 



o 



Here /, g are complex functions, u(t) denotes a sample path of the conventional 
Wiener process and dpo stands for the Wiener measure. Clearly, the kernel itself 
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can be explicitly written in terms of the conditional Wiener measure dp^^ pinned 
at space-time points (y, s) and (x, t), < s < t < T: 

k(y, s, x,t) = J exp[- j c(w(r), r)dr] dfi^(u) (4) 

As long as we do not impose any specific domain restrictions on the semigroup 
generator H(t), the whole real line R 1 is accessible to the process. Various choices 
of the Dirichlet |12|, §] boundary conditions can be accounted for by the formula (3). 
If we replace R 1 by any open subset Q C R 1 with the boundary dQ, it amounts to 
confining Wiener sample paths of relevance to reside in (be interior to) Q, which in 
turn needs an appropriate measure dp^ s l{u} G Q) in (4). This is usually implemented 
by means of stopping times for the Wiener process, |7|, 
Let f(x) and g(x) be two real functions such that: 

m T (x,y) = f(x)k(x,0,y,T)g(y) (5) 

defines a bi-variate density of the probability measure: 

myfi, B) = dx dym,T(x,y) (6) 
J a Jb 

i.e. a transition probability of the propagation from the Borel set A to the Borel set 
B to be accomplished in the time interval T. In particular, we need the marginal 
probability densities to be defined: 

p (x) = m T (x, n) , p T (y) = mr(fl, y) (7) 

where Vt C R is a spatial confining the process. 

Formulas (5), (6), can be viewed as special cases of (3), so establishing an appar- 
ent link between the Schrodinger problem and the Feynman-Kac kernels, together 
with the related parabolic equations (1). Assuming that marginal probability mea- 
sures (7) and their densities are given a priori, and a concrete Feynman-Kac kernel 
(4) (with or without Dirichlet domain restrictions) is specified, we are within the 
premises of the Schrodinger boundary data problem. 

Let Q = Q U dQ be a closed subset of R 1 , or R l itself. For all Borel sets (in 
the o"-field generated by all open subsets of Q) we assume to have known po(A) and 
Pt{B), hence the respective densities as well. If the integral kernel k(x,0,y,T) in 
the expression (5) is chosen to be continuous and strictly positive on Q, then the 
integral equations (7) can be solved || with respect to the unknown functions f(x) 
and g{y). The solution comprises two nonzero, locally integrable functions of the 
same sign, which are unique up to a multiplicative constant. 



i3,rara. 
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If, in addition, the kernel k(y,s,x,t), 0<s<t<Tisa fundamental solution 
16] of the parabolic system (1) on R 1 (i.e. is a function which solves the forward 



equation in (x,t) variables, while the backward one in (y,s)), then we have defined 
a solution of the system (1) by: 

u(x,t) = f(x,t) = / f(y)k(y,0,x,t)dy (8) 



v(x,t) = g(x,t) = J k(x,t,y,T)g(y)dy 

Moreover, p(x,t) = f(x,t)g(x,t) is propagated by the Markovian transition proba- 
bility density: 

p(y, s, x, t) = k(y, s, x, t) 9 Xj (9) 



p(x, t) = J p(y, s)p(y, s, x, t)dy 
< s < t < T 
d tP = A P - V{bp) 
Vg(x,t) 



b(x, t) 



g(x,t) 



the result, which covers all traditional Smoluchowski diffusions || 0. In that case, 
c(x, t) is regarded as time- independent, and the corresponding stochastic process is 
homogeneous in time. The Dirichlet boundary data can be implemented as well, thus 
leading to the Smoluchowski diffusion processes with natural boundaries, ||. Then, 
k(y, s, x, t) stands for an appropriate Green function of the parabolic boundary-data 
problem, with the property to vanish at the boundaries dQ of Q. 



Recently [^TJ , an extension of the above formalism was elaborated to encompass 
situations when the involved Feynman-Kac kernels are strictly positive and contin- 
uous, but not necessarily fundamental solutions of (1) nor even differentiable. They 
still give rise to (8) and (9) under suitable regularity conditions for solutions of the 
parabolic system (1). Their existence is not in conflict with the fact that k(y, s, x, t) 
itself needs not to be a solution of any differentiable equation. 

Let us also mention that for time-independent potentials, c(x, t) = c(x) for all 
t G [0,T], a number of generalizations is available |7|> §, |T§ [Lj, [Lj, |T8|-^T| to 
encompass the nodal sets of p(x) and hence of the associated functions f(x),g(x). 
The drift b{x) = X7lnp(x) = singularities do not prohibit the existence of 

a well defined Markov diffusion process (9), for which nodes are unattainable. In 
the considered framework they are allowed only at the boundaries of the connected 
spatial f2 confinig the process. 
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The problem of relaxing the strict positivity (and/or continuity) demand for 
Feynman-Kac kernels is nontrivial || |], [5] with respect to the eventual construction 
of the unique Markov process (9). To elucidate the nature of difficulties underlying 
this issue, let us consider quantally motivated examples of the parabolic dynamics 

(I)- 



2 Nonlinear parabolic dynamics with the funda- 
mental solution 

Let us choose the potential function c(x, t) as follows: 

c(x ' t) = 2(iT^-TT^ (10) 

for x G R 1 , t e [0, T\. In view of its local Holder continuity (cf. Ref. ||11|| ) with expo- 
nent one, and its quadratic boundedness, the fundamental solution of the parabolic 
system is known to exist |[16| , |22| , |23| , |24| . It is constructed via the parametrix method 
Among an infinity of regular solutions of (1) with the potential (10), we can in 
particular identify Jll] solutions of the Schrodinger boundary data problem for the 
familiar (quantal) evolution: 



p (x) = (27r)- 1/2 exp[-^] — > p(x,t) = [2rr(l + 1?)]-Wexp[-^^] (H) 



They read 

u{x, t) = f(x, t) = [2tt(1 + t 2 )V 1/A exp(-^Y^ + ^arctan t) (12) 

x 2 1 — t 1 

v(x, t) = g(x, t) = [2tv(1 + t 2 )Y l ^exp{— — - — ^ — -arctan t) 

and, while solving the nonlinear parabolic system (1) (with c = Ap 1 / 2 / p 1 ^ 2 ), in 
addition they imply the validity of the Fokker-Planck equation: 

p(x, t) = f(x, t)g(x, t) -> d t p = Ap - V(6p) 

b(x,t) = 2 =-—-x 13 

g(x,t) 1+t 2 

Notice that p(y, s, x, t) = k(y, s, x, t) is a fundamental solution of the first and 
second Kolmogorov (e.g. Fokker-Planck) equations in the present case. 
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Let us recall that a concrete parabolic system corresponding to solutions (12) 
looks badly nonlinear. Our procedure, of first considering the linear system (but with 
the potential "belonging" to another, nonlinear one), and next identifying solutions 
of interest by means of the Schrodinger boundary data problem, allows to bypass this 
inherent difficulty. In connection with the previously mentioned quantal motivation 
of ours, let us define g = exp(R + S),f — exp(R — S) where R(x, t), S(x, t) are real 
functions. We immediately realize that (1), (10) provide for a parabolic alternative 
to the familiar Schrodinger equation and its temporal adjoint: 

id t ip = -Aip (14) 

id t tp = Aip 

with the Madelung factorization ip = exp(R + iS), ip = exp(R — iS) involving the 
previously introduced real functions R and S. 



3 Nonlinear parabolic dynamics with unattain- 
able boundaries: the Green function 

Things seem to be fairly transparent when the parabolic system (1) allows for fun- 
damental solutions. However, even in this case complications arise if nodes of the 
probability density are admitted. The subsequent discussion has a quantal origin 
again, and comes from the free Schrodinger propagation (14) with the specific choice 
of the initial data: 

2 

M*) = * exp{~) — (15) 

2 

*p(x,t) = (2-r 1/4 irT ^ 7 , e ^[-i(TT^ ] 

such that our nonstationary dynamics example displays a stable node at x = for 
all times . 

The parabolic system (1) in this case involves the potential function: 

( «- A ^ 1/2 fo*) - x " 3 

WJ- p l,2 {X)t) - 2 (l+*2)2 l +t 2 ^> 



p(x, t) = (2tt)- 1/2 (1 + t 2 )- 3/2 x 2 exp[- 



x 2 



2(1 +t 2 ) J 

The polar (Madelung) factorization of Schrodinger wave functions implies: 

R(x,t) = lnp 1/2 (x,t) (17) 
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x > — > SYx, i) = (S+fx, t) = - — ~ — ~ — -arctan t 

4 1 + r 2 

t 3 

x < — > 5(x, i) = S-(x, t) = arctan t + n 

4 X ~\~ t 2 

Although S(x,t) is not defined at x — 0, we can introduce continuous functions 
/ = exp(R — S) and g = exp(R + 5") by employing the step function e(x) = if 
x > and e(x) = 1 if x < 0. Then, the candidates for solutions of the parabolic 
system (1) with the potential (16) would read: 

v(x, t) = g(x, t) = (27r) _1 / 4 (l + t 2 ) _3 / 4 |x| ex-p{—^—\ — \) exp[— -arctan t + 7re(x)] 



X 2 1 A- t 3 

w(x,t) =f(x,t) = (2tt)~ 1/4 (1 + t 2 )" 3/4 |x| exp( — ^^^) exp[-arctan t - vre(x)] 
For all x 7^ we can define the forward drift 



4 1 + t 2/ 2 
x 2 1 +f . r 3 

-JTTT 2)exp[ 2' 



6 x,t =2 , v y = - - x- 19 

v ' ' g(x,t) x 1 + t 2 y J 

which displays a singularity at x = 0. Nonetheless, (bp)(x,t) is a smooth function 
and the Fokker-Planck equation d t p = Ap — V(6p) holds true on the whole real 
line R 1 , for all t G [0, T]. Notice that there is no current through x = 0, since 
v(x,t) = 2VS(x,t) = vanishes at this point for all times. 

Our functions /(x, t),g(x, t) are continuous on R 1 , which however does not imply 
their differentiability. Indeed, they solve the parabolic system (1) with the potential 
(16) not on R l but on (— oo,0) U (0, +oo). Hence, almost everywhere on R 1 , with 
the exception of x = 0. 

An apparent obstacle arises because of this subtlety: these functions are not even 
weak solutions of (1), because of: 

/+oo f+OO 
d t f(x,t)(f)(x)dx+ / Vf{x,t)V<j)(x)dx+ (20) 

X f+OO 

- / c(x,t)f(x,t)(j>(x)dx 7^0 

Z J— oo 

for every test function such that 0(0) ^ 0, continuous and with support on a 
chosen compact set (e.g. vanishing beyond this set). 

One more obstacle arises, if we notice that c(x, t), (16) permits the existence 
of the unique, bounded and strictly positive fundamental solution for the parabolic 
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system (1). Then, while having singled out a fundamental solution and the bound- 
ary density data po(x), Pt(x) consistent with (16), we can address the Schrodinger 
boundary data problem associated with (2), (3): 

u(x,0) J k(x, 0,y,T)v(y,T)dy = p (x) (21) 

v(x,T) J k(y, 0,x,T)u(y,0)dy = p T (x) 

expecting that a unique solution u(x, 0),v(x,T) of this system of equations implies 
an identification u(x, 0) = f(x, 0) and v(x, T) = g(x, T). 

However, it is not the case and our f(x,t),g(x,t) do not come out as solutions 
of the Schrodinger problem, if considered on the whole real line R 1 , on which the 
fundamental solution sets rules of the game. 

Indeed, let us assume that (21) does hold true if we choose u(x, 0) = f(x, 0), v(x, T) = 
g(x,T), with / and g defined by (18). Since, in particular we have 

g(x, T) J k(y, 0, x, T)f(y, 0)dy = g(x, T)f(x, T) (22) 

then for there holds: 

f(x,T) = J k(y,0,x,T)f(y,0)dy (23) 

Both sides of the last identity represent continuous functions, hence the equality is 
valid point-wise (i.e. for every x). We know that f(y, 0) is continuous and bounded 
on R 1 , and k(y,0,x,T) is a fundamental solution of (1). Hence the right-hand-side 
of (23) represents a regular solution of the parabolic equation. Such solutions have 
continuous derivatives, while our left-hand-side function f(x,T) certainly does not 
share this property. Consequently, our assumption leads to a contradiction and (23) 
is invalid in our case. 

It means that the fundamental solution (e.g. the corresponding Feynman-Kac 
kernel) associated with (16) is inappropriate for the Schrodinger problem analysis, 
if the interpolating probability density is to have nodes (i.e. vanish at some points). 

In our case, x = is a stable node of p(x,t), and is a time-independent re- 
pulsive obstacle for the stochastic process. An apparent way out of the situation 
comes by considering two non-communicating processes, which are separated by the 
unattainable barrier at x — 0, [ |15| , 
The function 

f + (x, t) = (2tt)- 1/4 (1 + t 2 )" 3/4 x exp(-^-^-^r) exp{\arctan t) (24) 

4 1 + 1 2 



20, ft 25] 



8 



x G [0, oo) , t G [0, T] 



is a regular solution [TE\ of the first initial-boundary value problem for d t f = Af—cf 
specified by: 

f + (x,0) = (2n)- 1 / 4 xex P (-^) (25) 

/+(o,t) = o 

Then, instead of the fundamental solution, we need to utilize the Green function 
of the problem. To distinguish it from the fundamental Feynman-Kac kernel k we 
shall denote this Green function k + . Its existence is granted by the very existence 



of the fundamental solution for the considered potential (16), see Ref. fl38 [. 

The Green function k + (y, s,x,t) is a unique function such that for every con- 
tinuous on (0, oo) and with a compact support, the function: 

u(x,t) = k+(y,s,x,t)(j)(y)dy (26) 
J o 

is a solution of d t u = Au—c(x, t)u in (0, oo) x (s, T), with the properties: lim t i s u{x, t) 
4>(x) for all x G [0,oo), and u(0,t) = for all t G (s, T\. Moreover, for every 
(y, s) G (0, oo) x [0,T) the function k + is strictly positive in (0, oo) x (s, T) and 
k + (y,s,0,t) = for all t G {s,T). 

The uniqueness of solutions for the first initial-boundary value problem implies 
the validity of the semigroup composition rule for k + (Chapman-Kolmogorov iden- 
tity, which anticipates the Markov property of the constructed stochastic process). 

In this (uniqueness) connection a more detailed comment is necessary. 

For all s G [0, T), the function f + (x,s) can be uniformly approximated by a 
sequence of continuous functions <fi s n {x) such that, for each natural number n the 
support of (p s n is compact in (0, oo). There follows that the sequence of solutions of 
the first initial-boundary value problem d t u = Au — cu, given by: 

rod 

u n (x,t) = k + (y,s,x } t)(f) s n (y)dy (27) 
Jo 

is uniformly convergent to the solution f + (x,t). It implies that for any s < t we 
have _ 

f + (x,t)=l k+(y,s,x,t)f + (y,s)dy (2£ 

Now, let us consider 



o 



x 2 1 — t 3 

g + (x, t) = (2tt)~ 1/4 (1 + t 2 y 3/i x exp( — exp{--arctan t) (29) 

4 1 + t 2 2 
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which is the solution of the first initial-boundary value problem for the adjoint 
parabolic equation 

d t v = -Av + cv (30) 

g + (x,0) = (2ny 1 / 4 xexp(-^) 
g+(0,t) = 

for all t G [0, T\. Let k+ denotes the Green function of this adjoint equation. For 
every continuous function with a compact support in (0, oo) the formula 



v(y,s) = / kl.(x, t, y, s)(/)(x)dx (31) 
Jo 

with s < t, defines the solution of the first initial-boundary problem for the ad- 
joint equation. The previous arguments (at least for T < 1, modulo appropriate 
rescalings) apply in this well. We conclude that there holds 



9+(y, s )=l_ k* + (x,t,y,s)g + (x,t)dx (32) 

But, we have HT6|: 



kX(x,t,y,s) = k+(y,s,x,t) (33) 
for all x, y G (0, oo), and k+(x, t, 0, s) = 0. So, we can write: 

g+(y,s)= k+(y,s,x,t)g + (x,t)dx (34) 
Jo 

for y > 0, while g+(y, s) = if y = 0. 

All that finally allows us to introduce the transition probability density of the 
Markov process respecting the stable repulsive boundary at x = as follows: 

P+{y, s, x, t) = k + (y, s, x, t) 9+ Xj * (34) 

g+{y,-v 

y G (0, oo) , x G [0, oo) , < s < t < T 

The behaviour of p + at y = is to some extent irrelevant, and may involve a 
discontinuity. But, an innocent modification on the set of measure zero is allowed, 
and we choose p(0, s, x, t) = S(x). 

It completes the definition of the transition probability density of the Markov 
process, which is consistent with the dynamics of p(x, t). For all x G [0, oo), we have 
p{x,t) = j °° p+(y, s,x,t)p(y, s)dy, < s < t < T and also / °° p+(y, s,x,t)dx = 1 
for all y G [0, oo). 
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However, in view of b + (x, t) = 2 ^pfe- = ~ — x j^i which is singular at x = 0, 
the density p+(y, s,x,t) cannot by itself be a Green function for the associated 
Fokker-Planck equation d t p = Ap — V(bp), if considered on the whole of R 1 . The 
equation d t p + = Ap + — V(6p + ) holds true in the open set (0, oo) x (0,T). 

By combining the known results |L3|, [15], [I8|-||21|| about the unattainability of 
nodes by the diffusion process on R + (respectively, on we conclude that 

p+(y, s,x,t) respectively ) is a transition probability density of the diffusion 
with the density p(x,t), (16), forward drift b(x,t), (19) for which x = is an in- 
accessible repelling barrier. It remains in conformity with situations met in the 
conservative c(x, t) = c(x) cases, when an (ergodic p5| ) decomposition into the 
non-communicating due to nodes processes, is generic. 

On both semi-axes, the respective strictly positive, continuous (domain-restricted 
Feynman-Kac) kernels are defined almost everywhere, except for the barrier location 
x = 0, where they vanish. Thus, if considered on (— oo, 0] or [0, +oo), the respective 
integral kernels are non-negative , and no longer strictly positive. Moreover, they 
seem to need to be considered separately on R + and 

It is instructive to add that the existence of a node at time t = does not 
automatically imply its survival for times t > 0, and in reverse (while in the present 
context of non- negative kernels). 



4 The "Wiener exclusion" 

The conventional definition of the Feynman-Kac kernel (in the conservative case) 

exp[-t(-A + c)](y,x) = J exp[- c(u(T))dT]dp^(u) (35) 

comprises all sample paths of the Wiener process on R 1 , providing merely for their 
nontrivial redistribution by means of the Feynman-Kac weight exp[— f c(u(r))dr] 
assigned to each sample path u(s) : u(0) = y,u(t) = x. 

Assume that c(x) is bounded from below and locally (i.e. on compact sets) 
bounded from above. Then, the kernel is strictly positive and continuous |12| . 

For c = we deal with the conditional Wiener measure 

exp(tA)(y,x) = /$g[o;(s) € R 1 ; < s < t] = (36) 

fi[u(s) G R 1 ; w(0) = y, u(t) = x;0<s<t] 

pinned at space-time points (y, 0) and (x,t). 

The previous discussion indicates that i?_ is inaccessible for all sample paths 
originating from R + . In reverse, R + is inaccessible for those from R_. As well, we 
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may confine the process to an arbitrary closed subset fl C R 1 , or enforce it to avoid 
("Wiener exclusion" of Ref. |26|) certain areas in R l . 

In this context, it is instructive to know that [[T^] for an arbitrary open set fl, 
there holds: 

exp(tA n )(y,x) = /xg}[w(s) G fl, < s < t] (37) 

which is at the same time a definition of the operator — An, i.e. the Laplacian with 
Dirichlet boundary conditions, and that of the associated semigroup kernel. This 
formula provides us with the conditional Wiener measure which is confined to the 
interior of a given open set,p7L |28], B[. 



We can introduce an analogous measure, which is confined to the exterior of a 
given closed subset S C R . In case of not too bad sets (like an exterior of an 
interval in R l or a ball in R n , the corresponding integral kernel in known |12 to 
be positive and continuous. Technically, if S is a (regular) closed set such that the 
Lebesgue measure of dS is zero, then: 



exp(tA R \ s )(y, x) = fx%§[w(s) £S;0<s<t} 



(3* 



The Feynman-Kac spatial redistribution of Brownian paths can be extended to 
cases (37), (38) through the general formula valid for any f,g G L 2 (fl), where fl is 
any open set of interest (hence R\S, in particular): 



(f,exp(-tH n )g) = f(u(0))g(u(t))exp[- c(uj(T))dT}dfio(uj 



(39) 



It gives rise to the integral kernel comprising the restricted Wiener path integration, 
which is defined at least almost everywhere in x, y. Then, its continuity is not 
automatically granted. We can also utilize a concept of the first exit time Th for the 
sample path started inside fl (or outside S) 



T a (u)=inf[t > 0,X t (u) i ft] 
where X t is the random variable of the process. Then, we can write, K 

exp(-tHci)(y,x) = j exp[- ^ c(u(r))dT]diJ^[ijj)t < T n ] = 



(40) 



(41) 



exp[~ J Q c(w(r))tZr]d/ig; t) ; (w) 

It is an integration restricted to these Brownian paths, which while originating from 
y G fl at time t — are conditioned to reach x G fl at time t > without crossing 
(but possibly touching) the boundary dS of S. The contribution from paths which 
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would touch the boundary without crossing, for at least one instant s G [0, t] is of 
Wiener measure zero, ||27] . 



In case of processes with unattainable boundaries, with probability 1, there is 
no sample path which could possibly reach the barrier at any instant s < oo. 
The above discussion made an implicit use of the integrability property 

t 

c(uj(s))ds < oo (42) 

for uj G R l , < s < t., in which case the corresponding integral kernel (for bounded 
from below potentials) is strictly positive. Then, if certain inaccessible to 

the process, it occurs excusively [I%|]-[]2"I] due to the drift singularities, which are 



capable of "pushing" the sample paths away from the barriers. 

The previous procedure can be extended to the singular [p9|-|p7| potentials, 



which are allowed to diverge. Their study was in part motivated by the so called 
Klauder's phenomenon (and the related issue of the ground state degeneracy of 
quantal Hamiltonians), and had received a considerable attention in the literature. 

In principle, if S is a closed set in R 1 like before, and c(x) < oo for all x G 
Q = R\S, while c(x) = oo for x G S, then depending on how severe the singularity 
is, we can formulate a criterion to grant the exclusion of certain sample paths of 
the process and hence to limit an availability of certain spatial areas to the random 
motion. Namely, in case of (42) nothing specific happens, but if we have 

t 

c{uo(T))dr = oo (43) 

for lo(t) G S for some r G [0, t], then the "Wiener exclusion" certainly appears: we 
are left with contributions from these sample paths only for which (43) does not 
occur. Unless the respective set is of Wiener measure zero. 

The area Q comprising the relevant sample paths is then selected as follows: 

Q = b; f 'c(u(r))dT < oo] (44) 
Jo 

In particular, the criterion (44) excludes from considerations sample paths which 
cross 5* and so would establish a communication between the distinct connected 
components of Q. 

The singular set S can be chosen to be of Lebesgue measure zero and contain a 
finite set of points dividing R into a finite number of open connected components. 
With each open and connected subset Q C R 1 we can |29| associate a strictly positive 
Feynman-Kac kernel, which can be expected to display continuity. 
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Since the respective potentials diverge on S, their behaviour in a close neigh- 
bourhood of nodes is quite indicative. For, if us is a Wiener process sample path 
which is bound to cross a node at < s < t, then the corresponding contribution to 
the path integral vanishes. Such paths are thus excluded from consideration. If their 
subset is sizable (of nonzero Wiener measure), then the eliminated contribution 

I exp[- f c(uj s (T))dT}d^Uuj) = (45) 

is substantial in the general formula (35). 

At the same time, we get involved a nontrivial domain property of the semigroup 



generator H = — A + c resulting in the so called ground state degeneracy [29, 30, 31 



Let us recall (Theorem 25.15 in Ref. |[L2|| ) that if c is bounded from below and locally 
bounded from above, then the ground state function of H = — A + c is everywhere 
strictly positive and thus bounded away from zero on every compact set. 

5 Singular potentials and the ground state degen- 
eracy 

Our further discussion will concentrate mainly on singular perturbations of the har- 
monic potential. Therefore, some basic features of the respective parabolic problem 
are worth invoking. The eigenvalue problem (the temporally adjoint parabolic sys- 
tem now trivializes): 

-Ag + (x 2 — E)g = = Af — [x 2 - E)f (46) 

has well known solutions labeled by E n = 2n + 1 with n = 0, 1, 2, .... In particular, 
go(x) = fo(x) = -^j^exp{—^—) is the unique nondegenerate ground state solution. 
The corresponding Feynman-Kac kernel reads 

exp(-tH)(y,x) = k(y,0,x,t) = k t (y,x) = 

(n)-^(l - exp{-t)r^exp[- X ^- - ^"^ l (47) 

d t k = -A x k + (x 2 - l)k 

and the invariant probability density p(x) = f(x)g(x) = (7r) _1 / 2 exp(— x 2 ) is pre- 
served in the course of the time-homogeneous diffusion process with the transition 
probability density 

g (x^j 

p(y,s,x,t) = k t _ s (y,x)—- (48) 

9\y) 
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We have p(y, s, x, t) = p(y, 0, x, t — s). 

Notice the necessity of the eigenvalue correction of the potential both in (35) 
and (36), which is indispensable to reconcile the functional form of the forward drift 
b(x) = 2Vlng(x) = — 2x with the general expression for the corresponding (to the 
diffusion process) parabolic system potential 

1 b 2 

c = c ( X;t ) = d t lng + -(- + Vb) (49) 

which equals c(x) = x 2 — 1 in our case. 

Let us pass to the singular (degenerate) problems. 

Example 1: The potential 0: 

c(x) -x 2 + 72 - 7(7 + 1) - m 
C[X) ~ X + x 2 ^) \x\ 2 +~< |x|7 m 

with 7 > 0,x G R 1 , is singular at x = and is a well defined even function if 
otherwise. We can give [0] a solution to the stationary parabolic system 

-Ag + (c - l)g = = A/ - (c - 1)/ (51) 

in terms of 

g(x) = f(x)=exp[-(j^ + X ^)} (52) 

In this case, an invariant density up to normalization reads p(x) = (fg)(x) = g 2 (x) 
and is integrable on R 1 . It vanishes at x = and at both spatial infinities. 

By independent arguments ||15||-||21|| we know that a Markov diffusion process 



preserving p(x) can be consistently defined. The node x = is unattainable in view 
of the appropriate singularity of the forward drift: 

b(x) = 2Vln g{x) = sgn x - t — ~ %x (53) 



\x 



which pushes sample paths away from the node. Hence, there is no communication 
(realised by sample paths of the process) between R + and i?_ . Like in case of n > 1 
eigenfunctions of the harmonic oscillator, we deal with the totally disjoint (ergodic, 
|25|j ) components of the would-be global |2(J diffusion. 

This feature is nicely manifested in the apparent domain degeneracy of the as- 
sociated semigroup generator H = —A + (c — 1). Namely, Hg = is simply an 
eigenvalue problem. Let us define g + (x,t) = g(x,t) for x > and g+(x, t) = for 
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x < 0, and g-(x,t) = for x > while g_(x,t) = g{x,t) for x < 0. The same 
procedure can be repeated for / — > f±. 

The functions g + and g_ (/ + and /_ respectively) belong to L 2 (i? 1 ), and are 
orthogonal on R l , while corresponding to the same eigenvalue. It is an obvious 
spectral degeneracy of the respective generator H. As mentioned before, semigroups 
exp(-tH) with strictly positive kernels do not have |29[ generators with the ground 
state degeneracy. On the other hand, if 5* is the (singular) set of Lebesgue measure 



zero and R\S has m or more connected components, then there always exists |29 
a positive c(x) in Lj 0C (R\S) such that the ground state of H = —A + c is m-fold 
degenerate. 

This phenomenon we encounter in connection with (50). Recall that L\ oc com- 
prises equivalence classes of functions which are integrable on compact sets (e.g. 
locally). 

Example 2: The canonical (in the context of Refs. f2"5||-|p4j|) choice of the 
centrifugal potential: 

c E ( x ) = x 2 + ^-E (54) 

x A 



generates a well known spectral solution [33, 34 1 for Hg = [—A + cg(x)]g. The 
eigenvalues: 

£ n = 4n + 2 + (l + 8 7 ) 1/2 (55) 

with n = 0,1,2,... and 7 > -| (1 + 87) 1 / 2 > V%, are associated with the 
eigenfunctions of the form: 



2 



g n (x) = x^ +1 >/ 2 exp(~) L a n {x 2 ) (56) 



a = (1 + 87) 1/2 

n 

l* 2 ) = E 



n + a)\ (-x 2 Y 



u=0 



(n — is)\(a + v)\ v\ 



L%{x 2 ) = 1 , L^x 2 ) = -x 2 + a + l 

It demonstrates an apparent double degeneracy of both the ground state and of the 
whole eigenspace of the generator H. The singularity at x = does not prevent 
the definition of H = —A + x 2 + ^ since this operator is densely defined on an 
appropriate subspace of L 2 (R l ). This singularity is sufficiently severe to decouple 
(—00, 0) from (0, 00) so that L 2 (— 00, 0) and L 2 (0, 00) are the invariant subspaces of 
H with the resulting overall double degeneracy. 
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Potentials of the form, @ [26|, 



c{ x ) = x 2 + [dist(x, dtt)]- 3 (57) 

where dfl can be identified with dS, and S is a closed subset in R 1 of any (zero or 
nonzero) Lebesgue measure, have properties generic to the Klauder's phenomenon. 
Because the Wiener paths are known to be Holder continuous of any order | — e, e > 
and of order | in particular, there holds Jq c{u{r))dr = oo if lu(t) E S for some r. 
Conversely, J c (u;(r))(ir < oo if u never hits S 1 . This implies that the relevant 
contributions to: 

(/, exp[-t(-A + c)]g) = J J(u(0))g(u(t))exp[- jf* c(w(r))dr]^ (w) (58) 
come only from the subset of paths defined by 

Q t = [u- f c{u{r))dr < oo] (59) 
jo 

The above argument might seem inapplicable to the centrifugal problem. However 
it is not so. In the discussion of the divergence of certain integrals of the Wiener 
process, in the context of Klauder's phenomenon, it has been proven |56] that for 
almost every path from x — 1 to x = — 1 (crossing the singularity point x = 0) there 
holds f£*g |c<j(r)| _1 (ir = oo for any 5 > 0. 

To be more explicit: if T\ = t\{oj) is the first time such that the Wiener process 
W(t) = W(t,u) attains the level (location on R 1 ) W{t\) = 1, then the integral over 
any right-hand-side neighbourhood (ji, T\ + 5) of T\ diverges: 

[ T1+S c(uj(t) - l)dt = oo (60) 

if J^i c(x)dx = oo. 

In case of the left-hand-neighbourhood of T\, we have 

f T1 c{u(t) - l)dt = oo (61) 

if xc(x)dx = oo. 

All that holds true in case of the centrifugal potential, thus proving that the 
only subset of sample paths, which matters in (58) is (59). Obviously, Q t does not 
include neither paths crossing x = nor those which might hit (touch) x = at 
any instant. The singularity is sufficiently severe to create an unattainable repulsive 
boundary for all possible processes, which we can associate with the spectral solution 
(55),(56). 
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6 The singular potential in action 



After the previous analysis one might be left with an impression that the appearence 
of the stable barrier at x = persisting for all t G [0,T], is a consequence of the 
initial data choice ^o(O) = for the involved quantum Schrodinger picture dynamics. 
In general it is not so. For example, ipo(x) = x 2 exp(— x 2 /4) which vanishes at x = 0, 
does not vanish anymore for times t > of the free evolution. 

On the other hand, somewhat surprisingly from the parabolic (intuition) view- 
point, the node can be dynamically developed from the nonvanishing initial data 
and lead to the nonvanishing terminal data. 

Let us consider a complex function: 

*(,,«) = (1 + !i )-v^ p| __£l_ | ,_£!_ + (62) 

2 2 

which solves the free Schrodinger equation with the initial data ip(x, 0) = ^exp{— ^-J. 
It vanishes at x = exclusively at the initial instant t = of the evolution. 
Obviously, there is nothing to prevent us from considering 

ift(x,t)=il)(x,t-a) (63) 

for a > 0. It solves the same free equation, but with nonvanishing initial data. How- 
ever, the node is developed in the course of this evolution at time t = a and instanta- 
neously desintegrated for times t > a. Here, the Schrodinger boundary data problem 
would obviously involve two strictly positive probability densities po(x) = \^(x, 0)| 2 
and Pt(x) = \^(x,T)\ 2 , T > a. It would suggest to utilize the theory fllTf , based on 
strictly positive Feynman-Kac kernels, to analyze the corresponding interpolating 
process. However, this tool is certainly inappropriate and cannot reproduce the a 
priori known dynamics, with the node arising at the intermediate time instant. 

To handle the issue by means of a parabolic system, which we can always asso- 
ciate with a quantum Schrodinger picture dynamics, let us evaluate the potential 
c(x,t) appropriate for (1). 

In view of 

2 4 

p(x, t) = const (1 + t 2 )- 5 / 2 exp[- X + ] - x 2 t 2 + t 2 (l+ t 2 )] (64) 
we have (while setting w 1 ^ 2 (x, t) — — x 2 t 2 + t 2 {\ + t 2 )}): 
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1 x 2 1 3x 2 -2t 2 x 1 

4(l+t 2 ) 2 ~ 2f- t 2 X 2 +t 2 {l + t 2 ) ~ 2(1 + t 2 ) + 

1 3x 2 -It 2 1 , x 3 - 2t 2 x 



2 f - t 2 X 2 + t 2 (l + t 2 ) 4 V ^ - £2 X 2 + f 2(! + t 2) ' 

The expression looks desparately discouraging, but its t J, (i.e. the initial data ) 
limit is quite familiar and displays a centrifugal singularity at x = 0: 

, N Ap x / 2 (x,0) x 2 2 5 

Since the original, dimensional expression for the centrifugal eigenvalue problem 
is §3|: 



("^A + ^x 2 + l)g = Eg (67) 

S n = m[2n + l + i(l + 8 7 ) 1 / 2 ] 

with n = 0,1,..., an obvious adjustment of constants m = 1/2,7 = 1 allows to 
identify E = 5/2 as the n = eigenvalue of the centrifugal Hamiltonian H = 

"A + ^ + J, 

A peculiarity of the considered example is that it enables us to achieve an ex- 
plicit insight into an emergence of the centrifugal singularity and its subsequent 
destruction (decay) for times t > a, due to the free quantum evolution. 

2 2 

In view of the degeneracy of the ground-state eigenfunction \exp{—^-) of the 
centrifugal Hamiltonian, we deal here with the gradually decreasing communication 
between R + and _R_, which results in the emergence of the completely separated 
(disjoint) sets (— oo, 0) and (0, +oo) at t = a, followed by the gradual increase of 
the communiaction for times t > a. By "communication" we understand that the 
set of sample paths crossing x = forms a subset of nonzero Wiener measure. 

It also involves a generalisation (cf. also Refs. 0, |37], 38| ) to time- dependent 
Feynman-Kac kernels: 

(/, exp[- f s H(r)dr}g) = J f(u(s))g(u(t))exp[- jT c(uj(t), r)dr] d^{oo) (68) 

Qs,t=[w, [ c(u (t), t) dr < oo] 



The finiteness condition f 8 c(uj(t), r)dr < oo, surely does not hold true, [Q, if 5 > 
is sufficiently small, cf. (60), (61). 
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Let us mention that a few interesting mathematical questions are left aside in 
the present paper. They deserve a separate, thorough study on their own. For 
example, even in case of conventional Feynman-Kac kernels, the weakest possible 
criterions allowing for their continuity in spatial variables are not yet established. 
An issue of the continuity of the kernel in case of general singular potentials, needs 
an investigation as well. 



7 Nonnegative solutions of parabolic equations and 
the Schrodinger boundary data problem accord- 
ing to R. Fortet 

As emphasized before, one of motivations for our analysis was the quantal observa- 
tion (exploiting the Born statistical interpretation postulate as the principal building 
block of the theory) that the temporally adjoint pair of Schrodinger equations: 

idtif>{x, t) = -DAip(x, t) + —^V(x)if>{x, t) (69) 

id$(x,t) = DA$(x,t) - ——V(xj$(x,t) 

2mD 

by means of the polar decomposition ip = exp(R + iS),ip = exp(R — iS) can be 
transformed into the (hopelessly looking at the first glance) nonlinearly coupled 
parabolic system of the form (1): 

OA = DA9* - ^(2Q - m (70) 
d t e = -DA9 + -^—(2Q-V)9 



Q = 2mD 



Ap 1 ' 2 

,1/2 



p(x,t) = 9*(x,t)9(x,t) = ij){x,t)il){x,t) 

for real functions 9 = exp(R + S),9* = exp(R — S). In the above h/2m = D can be 
set to restore the traditional notation. 

While searching for a probabilistic meaning of the system (70), we had in fact 
assumed (see also Refs. [Q, [11]) to have in hands a solution of (69), so that p(x,t) 
was known. Our next step amounts to replacing the nonlinear parabolic system (70) 
by a linear one (1), (with D — m — 1) where for each given functional choice of 
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c(x, t) = 2^j)i^Q — V), all allowed solutions u(x, t),v(x, t) were sought for, including 
the fundamental one. At this point the crucial role of the respective Feynman-Kac 
kernel was disclosed. Effectively, we have invoked the Schrodinger boundary data 
problem to pick up the unique solution from among an infinity of others, such that 
the arising probability measure dynamics (if any) is consistent with the a priori 
known probability density boundary data. 

In this indirect way, a definite justification for the probabilistic significance of 
the nonlinear parabolic system has been achieved. The procedure proves to be 
consistent, thus allowing to investigate fundamental solutions, Green functions and 
other solutions of (70). Even, if viewed as the nonlinearly coupled parabolic system 
per se. See e.g. also Refs. |], [39| for a discussion of the "propagation of chaos" in a 
system of interacting (coupled) diffusion processes 

If the Feynman-Kac kernels are strictly positive, the respective solutions are 
strictly positive as well, except for the boundaries dfl of the spatial area confining 
the process. Once we admit the nonnegative Feynman-Kac kernels, we fall into an- 
other theoretical framework, this of nonnegative solutions of linear (and nonlinear) 
parabolic equations, PD[ . Since we address the general time-independent potentials, 
and as we have seen the time-dependent domain properties of semigroup generators 
are involved, we touch upon an almost undeveloped theory of construction of Markov 
processes associated with time-dependent Dirichlet forms and spaces, ]|l|, |42| . The 



standard theory of forms pi] does not work in case of time-inhomogeneous evolu- 
tions. 

On the other hand, the Schrodinger problem itself needs an extension to non- 
negative Feynman-Kac kernels. The strongest (uniqueness of solution) result was 
established || for strictly positive kernel functions, a demand which needs to be 
relaxed for our purposes. 

Hence, our major assumption must be that the kernel is a nonnegative and 
continuous function. One can try to relax the continuity condition as well. Moreover, 
it appears that the kernel may not be a function, and nonetheless one can expect 
that the main features of our analysis would persist. 

Beurling |4j] has attempted to relax the strict positivity condition in more than 
one spatial dimension, with a partial success only. An earlier analysis due to Fortet 
|| is of particular importance in our one- dimensional context. He addressed an 
issue of the existence and uniqueness of nonnegative solutions of the Schrodinger 
boundary data problem, under an assumption that the kernel is continuous and 
nonnegative in one spatial dimension. 

Fortet 's integral equations: 

Pi{x) = f{x) k(x,y)g(y)dy (71) 
Jn 2 
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P2{x)=g(y) f(x)k(x,y)dx 

J f2i 



where Sli, Sl 2 are finite (or not) intervals in R 1 , and J Ui pi(x)dx = Jf l2 P2(y)dy > 0, 
while k(x,y) > 0, p\(x) > 0,p 2 (y) > 0, were to be solved with respect to the 
unknown functions f(x),g(y), defined respectively on Q 1 ,Q 2 . 

All functions k(x,y), pi(x), p 2 (y) are by assumption real and measurable (and 
integrable) for x e Qi, y e Q2- There are however the additional assumptions which 
must be respected: 

(i) k(x,y) is continuous, bounded from the above and nonnegative almost every- 
where in Qi x Q 2 , i-e. except for a set S of measure zero comprising both x's and 
y's in R 1 , 

(ii) pi(x) and p 2 (y) are continuous, 

(iii) let A be a closed interval in Q 1 . For a nonnegative continuous in Q 2 func- 
tions g(y) we demand that if the integral G(x) = J n2 k(x,y)g(y)dy is finite almost 
everywhere on an open subset of Qi containing A, then this integral is uniformly 
convergent on A. Analogously with respect to y G B C Q, with f(x) — > -F(y) on fi 2 . 

Under these hypotheses, the integral equations (71) admit a unique solution 
given in terms of two functions: 

(1) f(x) which is strictly positive and continuous almost everywhere, except for the 
set of zeros of pi(x), 

(2) g(y) which has the same zeros as f(x), is strictly positive almost everywhere 
and measurable. 

Hence, g(y) is not necessarily continuous, and one should realise that we must 
have granted the existence of Vln g(x,t) as the drift field. If the function g(x) 
is continuous, all that fits to our previous discussion. However, even in this case 
some additional restrictions on k(x, y) are necessary to guarrantee a differentiability 
of f(x,t) = J k(y,0,x,t)f(y)dy and g(x,t) = J k(x,t,y,T)g(y)dy, and make them 
solutions of the time-adjoint parabolic system once we set k(x, y) = k(y, s,x,t),0 < 
s < t < T and select the appropriate Feynman-Kac kernel. 

In connection with the previous centrifugal example, the above conditions appear 
to be too restrictive. Let us therefore invoke another result due to Fortet. Namely, 
if the above condition (iii) is replaced by the demand: 



/ 



P2{V) dy < 00 (72) 



^2 [/m k{z,y)p l {z)dz 

then a unique nonnegative solution of the integral equations (71) comprises a con- 
tinuous function f(x) whose zeros coincide with those of pi(x). The function g(y) is 
measurable and has zeros of p 2 (x), which are not necessarily in common with those 
of pi. 
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This result opens a number of interesting propagation scenarios, and deserves a 
careful analysis (with the generalization prospects) in higher dimensions, cf. also for 
a related discussion in Refs. [H3L [PJ. Notice, that our centrifugal example exhibits 



in its simplest version an intriguing feature of Fortet's analysis: we are capable 
of producing probability densities p\{x) (initial) and P2(y) (terminal), which have 
non- coinciding sets of zeros. The inquiry into the corresponding Schrodinger's in- 
terpolating dynamics is quite an appealing problem. 
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